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Abstract 

We consider possible signatures for the so called hybrid gravity within the Galactic Central Parsec. This modified 
theory of gravity consists of a superposition of the metric Einstein-Hilbert Lagrangian with an f(R) term constructed 
a la Palatini and can be easily reduced to an equivalent scalar-tensor theory. Such an approach is introduced in order 
to cure the shortcomings related to f(R) gravity, in general formulated either in metric or in metric-affine frameworks. 
Hybrid gravity allows to disentangle the further gravitational degrees of freedom with respect to those of standard 
General Relativity. The present analysis is based on the S2 star orbital precession around the massive compact dark 
object at the Galactic Centre where the simulated orbits in hybrid modified gravity are compared with astronomical 
observations. These simulations result with constraints on the range of hybrid gravity interaction parameter (pQ, show¬ 
ing that in the case of S2 star it is between -0.0009 and -0.0002. At the same time, we are also able to obtain the 
constraints on the effective mass parameter and found that it is between -0.0034 and -0.0025 AU“’ for S2 star. 
Furthermore, the hybrid gravity potential induces precession of S2 star orbit in the same direction as General Rela¬ 
tivity. In previous papers, we considered other types of extended gravities, like metric power law f{R) oc R" gravity, 
inducing Yukawa and Sanders-like gravitational potentials, but it seems that hybrid gravity is the best among these 
models to explain different gravitational phenomena at different astronomical scales. 
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1. Introduction 

The existence of different anomalous astrophysical 
and cosmological phenomena like the cosmic accel¬ 
eration, the dynamics of galaxies and gas in clus¬ 
ters of galaxies, the galactic rotation curves, etc. re¬ 
cently boosted the growth of several long-range modi¬ 
fications of the usual laws of gravitation. These men¬ 
tioned phenomena did not find satisfactory explana¬ 
tions in terms of the standard Newton-Einstein grav¬ 
itational physics, unless exotic and still undetected 
forms of matter-energy are postulated: dark matter and 
dark energy. A recent approach is to try to explain 
these phenomena without using new material ingredi¬ 
ents like dark matter and dark energy, but using well- 
motivated generalization and extensions of General Rel¬ 
ativity (GR). Several alternative gravity theories have 
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Newtonian gravity in the framework of the weak field 
limit ll26l] of fourth-order gravity theory have been pro¬ 
posed and constraints on these theories have been dis¬ 
cussed ES [III [11 Sill 111 [si 111 ESS. 

The philosophy is to search for alternative form of 
gravity, i.e. of the Einstein-Hilbert theory, so that 
such modifications could naturally explain some astro- 
physical and cosmological phenomena without invok¬ 
ing the presence of new material ingredients that, at the 
present state of the art, seem hard to be detected. Be¬ 
sides, this extended approach can be connected to effec¬ 
tive theories that emerge both from the quantization on 
curved spacetimes and from several unification schemes 

fiSli. 
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The simplest extension of the Einstein-Hilbert action 
is based on straightforward generalizations of the Ein¬ 
stein theory where the gravitational action (the Einstein- 
Hilbert action) is assumed to be linear in the Ricci cur¬ 
vature scalar R. If this action consists in modifying the 
geometric part considering a generic function f{R), we 
get so called f{R) gravity which was firstly proposed 
in 1970 by Buchdahl |39]. Generally, the most seri¬ 
ous problem of f(R) theories is that these theories can¬ 
not easily pass the standard Solar System tests SEil. 
However, there exists some classes of them that can 
solve this problem It^] . It can be shown that f(R) the¬ 
ories, in principle, could explain the evolution of the 
Universe, from a matter dominated early epoch up to 
the present, late-time self accelerating phase. Several 
debates are open in this perspective ll43ll44l 4^ but 
the crucial point is that suitable self-consistent model 
can be achieved. f{R) theories have also been studied in 
the Palatini approach, where the metric and the connec¬ 
tion are regarded as independent fields 14711 . Metric and 
Palatini approaches are certainly equivalent in the con¬ 
text of GR, i.e., in the case of the linear Einstein-Hilbert 
action. This is not so for extended gravities. The Pala¬ 
tini variational approach leads to second order differen¬ 
tial field equations, while the resulting field equations 
in the metric approach are fourth order coupled differ¬ 
ential equations. These differences also extend to the 
observational aspects. 

A novel approach, that consists of adding to the 
metric Einstein-Hilbert Lagrangian an f{R) term con¬ 
structed within the framework of the Palatini formal¬ 
ism, was recently proposed |4^ 4^ The aim of 
this formulation is twofold; from one side, one wants to 
describe the extra gravitational budget in metric-affine 
formalism, from the other side, one wants to cure the 
shortcomings emerging in f(R) gravity both in metric 
and Palatini formulations. In particular, hybrid gravity 
allows to disentangle the metric and the geodesic struc¬ 
tures pointing out that further degrees of freedom com¬ 
ing from f(R) can be recast as an auxiliary scalar field. 
In such a case, problems related to the Brans-Dicke-like 
representation of f(R) gravity in terms of scalar-tensor 
theory (the so called O’Hanlon transformation) are im¬ 
mediately avoided (see if^ for details and the discus¬ 
sion in Sec. 2). Due to this feature, the scalar-tensor 
representation of hybrid gravity results preferable with 
respect to other scalar-tensor representations of gravi¬ 
tational interaction. As byproducts, the appearance of 
ghosts is avoided and the correct weak field limit of f(R) 
gravity with respect to GR is recovered. Eurthermore, 
several issues related to the galactic dynamics, the for¬ 
mulation of the virial theorem in alternative gravity, the 


dark energy behavior seem to be better addressed than 
in f(R) gravity considered in both metric or Palatini for¬ 
mulations. In summary, the hybrid metric-Palatini the¬ 
ory opens up new possibilities to approach, in the same 
theoretical framework, the problems of both dark en¬ 
ergy and dark matter disentangling the extra degrees of 
freedom of gravitational field with respect to GR. Eor 
a brief review on the hybrid metric-Palatini theory, we 
refer the reader to iH. 

In this perspective, star dynamics around the Galactic 
Centre could be a useful test bed to probe the effective 
gravitational potentials coming from the theory. In par¬ 
ticular, S-stars are the young bright stars which move 
around the centre of our Galaxy ^ ^ .53. .571 

where the compact radio source Sgr A* is located. Eor 
more details about S2 star see references OSSIdTII . There 
are some observational indications that the orbits of 
some of them, like S2, could deviate from the Keple- 


rian case 11-54 l59l] . but the current astrometric limit is 
not sufficient to unambiguously confirm such a claim 

iSS. 

Here we study a possible application of hybrid mod¬ 
ified gravity within Galactic Central Parsec, in order 
to explain the observed precession of orbits of S-stars. 
This paper is a continuation of previous studies where 
we considered different extended gravities, such as 
power law f(R) gravity ms, f(R, (p) gravity imply¬ 
ing Yukawa and Sanders-like gravitational potentials in 
the weak field limit Results obtained using 

hybrid gravity point out that, very likely, such a the¬ 
ory is the best candidate among those considered to ex¬ 
plain (within the same picture) different gravitational 
phenomena at different astronomical scales. 

More details about hybrid gravity you can find in 
|47, ^ M]. It is shown in M] that this type 

of modified gravity is coherently addressing the Solar 
System issues, and motivations for addressing them are 
discussed in details in M- 

The modified theory of gravity needs to be con¬ 
strained at different scales; at laboratory distances, at 
Solar system, at galaxies, at galactic clusters and at cos¬ 
mological scales. Obtaining constraints at any of these 
scales is a fundamental issue to select or rule out mod¬ 
els. In particular, it is important to investigate gravity 
in the vicinity of very massive compact objects because 
the environment around these objects is drastically dif¬ 
ferent from that in the Solar System framework. The S2 
star orbit is a unique opportunity to test gravity at the 
sub-parsec scale of a few thousand AU. Eor example, 
gravity is relatively well constrained at short ranges (es¬ 
pecially at sub-mm scale) by experimental tests, how¬ 
ever for long ranges further tests are still needed (see 
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Figures 9 and 10 from fill] for different ranges). It is 
worth stressing that a phenomenological approach can 
be useful in this context. In particular, the motion of 
S2-star is a suitable tool to test alternative gravity. For 
the reasons that we will discuss in detail below, hybrid 
gravity is a reliable paradigm to describe gravitational 
interaction without considering dark energy and dark 
matter. Specifically, the massive compact object inside 
the Galactic Center is surrounded by a matter distribu¬ 
tion and deviations of S2-star motion from the Keple- 
rian orbit are observed in detail. These deviations can 
be triggered both by the masses of the surrounding bod¬ 
ies and by the strong field regime at the Galactic Cen¬ 
ter. This peculiar situation constitutes a formidable op¬ 
portunity to test theories of gravity. However, it is im¬ 
portant to stress that numerical results reported here by 
comparing models with astronomical observations, rep¬ 
resent only upper bounds for the precession angle on the 
deviation from GR. More accurate studies will be neces¬ 
sary in future work to better constrain dynamics around 
the Galactic Center. 

The present paper is organized as follows: in Sec. 2 
we sketch the theory of hybrid gravity. In Sec. 3 we 
describe our simulations of stellar orbits in the gravita¬ 
tional potential derived in the weak field limit of hybrid 
gravity and we describe the fitting procedure to simu¬ 
late orbits with respect to astrometric observations of 
S2 star. Results are presented in Sec. 4. Conclusions 
are drawn in Sec. 5. 


2. Hibrid metric-Palatini gravity 


In this Section, we present the basic formalism for 
the hybrid metric-Palatini gravitational theory within 
the equivalent scalar-tensor representation (we refer the 

reader to if^ EE E2. E3] for more details). The f(R) 

theories are the special limits of the one-parameter class 
of theories where the scalar field depends solely on the 
stress energy trace T (Palatini version) or solely on the 
Ricci curvature R (metric version). Here, we consider a 
one-parameter class of scalar-tensor theories where the 
scalar field is given as an algebraic function of the trace 
of the matter fields and the scalar curvature 11621] : 


f(R) gravity in the Palatini and the metric formalism, re¬ 
spectively. For any finite value of Q. 4 , its value depends 
both on matter and curvature. In the limit Qyi —> cx) 
the propagating mode is given solely by the curvature, 
0(R, T) —> <p(R), and in the limit —> 0 solely the 

matter fields 0(R, T) 4>{T). In the general case, the 
field equations are fourth order both in the matter and in 
the metric derivatives as we will show below. 

More specifically, the intermediate theory with = 
1 and D = 4, corresponds to the hybrid metric-Palatini 
gravity theory proposed in where the action is 

given by 


5 = J [r + /(R) + 2k^J:„] . (2) 

where = SttG, R is the Einstein-Hilbert term, R. = 
g^''Rfiv is the Palatini curvature with the independent 
connection F^j^ as 




(3) 


The Palatini-Ricci tensor is 


R = f" - f“ + N _ f“ ('41 

— A Q, 1 y ^ 1 1 1 Q,y • 


fia,v 


^ av • 


Varying the action given with respect to the metric, one 
obtains the field equations 


G^y + F(R)R^, - ^f(R)g^y = , (5) 

where the matter stress-energy tensor is 


2 6 ( 


( 6 ) 


The independent connection is compatible with the met¬ 
ric F{R)g^y, conformal to g^y, with the conformal factor 
given by F{R) = df{R)ldR. This fact gives 


R^v = V + |p4jF(R),^F(R),y 

- Fk>'^,F('R).y - 5^^^vV„V“F(R) . 


(7) 



D-l 

2 (D- 2 )(Qa- 0 ) 


- V(<f>) 


( 1 ) 

The theories can be parameterized by the constant Q^. 
The limiting values = 0 and correspond to 

scalar-tensor theories with the Brans-Dicke parameter 
0 ) - -{D- l)l(D-2) and w = 0. These limits reduce to 


The Palatini curvature R is obtained from the trace of 
the field equations (|5]l, which is 

F(R)R - 2f{R) +R = X. ( 8 ) 

R can be algebraically expressed in terms of X if f{R) 
is analytic. In other words, the variable X measures how 
the theory deviates from GR trace equation R - -k^T. 
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We can express the field equations Q in terms of the 
metric and X as 


G 


fiv 


+ 


+ 


^f(X)g^y - F(X)R^y + F'(XW^X,y 
i [f'(X)V„V“X + F"(X) mf] g^v 


F”{X) 


3 {F'{X)f 
2 F{X) 


X,^X,y+K%ym 


being (dX)^ - The trace of the field equations is 

now 


T'(X)V„V“X + 


F”(X) ■ 


1 (F\X)Y 


2 F(X) 


mf 


1 


+ -[X + 2f{X)-F{X)R\^0, (10) 


while the relation between the metric scalar curvature R 
and the Palatini scalar curvature R. is 


RiX) = R 


F'jxW , V„V°'F(X) 
F{X) I F(X) 


( 11 ) 


which can be obtained by contracting Eq. 0 . As for 
pure metric and Palatini cases |01, the action 0 for the 
hybrid metric-Palatini theory can be recast into a scalar- 
tensor theory by an auxiliary field A such that 


5 = d^xY^[R + f(A) + fA(R-A)] + S„,, 

( 12 ) 

where /a = df /dA and Sm is the matter action. Rear¬ 
ranging the terms and defining (p = /a, yi<t>) - A/a - 
f{A), Eq. (fT^ becomes 


5 = ^ J" d'^x [R + cpR- V(p)] +Sm ■ 


(13) 


The variation of this action with respect to the metric, 
the scalar p and the connection leads to the field equa¬ 
tions 

R^y + (l)'R^,y-^-{R+<pR-V)g^y = ^%v,(14) 

R-V^ = 0, (15) 

= 0 , (16) 

respectively. The solution of Eq. (fTfil l implies that the 
independent connection is the Levi-Civita connection of 
a metric = <Ps^lV 7 that is we are dealing with a bi¬ 
metric theory and R^y and are related by 


R,iv = Riav + ^df,<pdyp - 7 ( V;jVy(^ + ^g^yVaX^fp 


(17) 


which can be used in the action (flTt to obtain the fol¬ 
lowing scalar-tensor representation 


S = 




(l+cp)R+—d^,pdf^,p-V{,p) 
Zcp 


+S n 
(18) 


We have to stress that, by the substitution <p —> 
-{K(pYl6, the action (fTSl) reduces to the case of a confor¬ 
mally coupled scalar field with a self-interaction poten¬ 
tial. This redefinition makes the kinetic term in the ac¬ 
tion (fTsT i the standard one, and the action itself becomes 
that of a massive scalar-field conformally coupled to the 
Einstein gravity. Of course, it is not the Brans-Dicke 
gravity where the scalar field is massless. 

As discussed above, in the limit Q.a —> 0, the the¬ 
ory (fTSl l becomes the Palatini-ZCR) gravity, and in the 
limit Q/i ^ oo it is the metric f{R) gravity. Apart from 
these singular cases, any theory with a finite Q .4 is in 
the ’’hybrid” regime, which from this point of view pro¬ 
vides a unique interpolation between the two a priori 
completely distinct classes of gravity theories. 

Using Eq. (fTTl) and Eq. (fTSl) in Eq. (fT4li . the metric 
field equations are 


(1 -H (p)Rf,y = lY {r^y - ig^yT j 

+-gnv{y + v„v“0) 

+yY^v(p-^df,(pdy(p, (19) 
l(p 

and then the spacetime curvature is sourced by both 
matter and scalar field. The scalar field equation can 
be manipulated in two different ways that illustrate how 
this theory is related with the w = 0 and w - -3/2 
cases, which corresponds to the metric and Palatini 
scalar-tensor representations of /(R)-gravity | 3 l respec¬ 
tively. Considering the trace of Eq. dll with g^'’, we 
find -R - (pR + 2V - k^T, and using Eq. (fT3T l. it is 


2V-(py^ = k^T+R . (20) 


Similarly as in the Palatini case (w = -3/2), this equa¬ 
tion says that the field (p can be expressed as an algebraic 
function of the scalar X = k^T + R, i.e., (p - (p{X). In 
the pure Palatini case, however, p is just a function of T. 
Therefore the right-hand side of Eq. (US contains mat¬ 
ter terms associated with the trace T, its derivatives, and 
also the curvature R and its derivatives. In other words, 
this theory can be seen as a higher-derivative theory in 
both matter and metric fields. However, such an inter¬ 
pretation can be avoided if R is replaced in Eq. (l20l l with 
the relation 
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( 21 ) 















together with K - V^. One then finds that the scalar 
field dynamics is given by a second-order equation that 
becomes, for Qa = 1 , 




0[2y-(i 


“ 3 ’ 

( 22 ) 

which is a Klein-Gordon equation. This result shows 
that, unlike in the Palatini case (w = -3/2), the scalar 
field is dynamical. In this sense, the theory is not af¬ 
fected by the microscopic instabilities that arise in Pala¬ 
tini models (see ll47n for details). 


3. The weak field limit and the fitting procedure 


In the weak field limit and far from the sources, the 
scalar field behaves as (p(r) ^ (po + (2G(f>oM 
the effective mass is defined as 


ml = (2y - - 0(1 + 0)V00)/3|^^^^ , (23) 


where 0 o is the amplitude of the background value of 0 . 
Furthermore V, and are respectively the potential 
and its first and the second derivatives with respect to 0 . 
The metric perturbations yield 




2GeffM To 

, 

r 1 -t- 00 6 
/ 2yGeffM yp 

\ r 1 +00 6 / '' 


(24) 


where yo is the minimum of the potential V. The effec¬ 
tive Newton constant Geff and the post-Newtonian pa¬ 
rameter y are defined as 

Geff = -^[1-(0o/3)e-'”^''] , 

1 +00 

^ l+(0o/3)c-^'' 


The coupling of the scalar field to the local system de¬ 
pends on 00 . If 00 ^ 1, then Geff ~ G and y » 1 regard¬ 
less of the value of This is in contrast with the result 
obtained in the metric version of f{R) theories. For suf¬ 
ficiently small 00 , this modified theory allows to pass 
the Solar System tests, even if the scalar field is very 
light 05111 . According to these considerations, the lead¬ 
ing parameters are and 0o. Their value give both an 
estimation of the deviation with respect to GR and how 
the affine contribution (i.e. the Palatini term) is relevant 
with respect to the metric f{R) gravity. Constraining 
both of them by observations gives immediately infor¬ 
mation on the hybrid gravity. Starting from the above 


results, the modified gravitational potential can be writ¬ 
ten in the form: 


0(0 = [1 - (0o/3)e-'”^''] Mir. (26) 

1 +00 

An important remark is necessary at this point. We have 
not chosen the form of y(0) since the only requirement 
is that the scalar field potential is an analytic function of 
0. In such a case, the effective mass (|2^ can be always 
defined. Clearly, the aim is to derive specific forms of 
the potential starting from the observations. This means 
a sort of ’’inverse scattering procedure” by which the 
y(0) potential can be reconstructed from the observed 
values of the parameters M, 0o, and y. 

To this end, let us use eq. (l26l l to simulate orbits 
of S2 star in the hybrid modified gravity potential and 
then we compare the obtained results with the set of S2 
star observations obtained by the New Technology Tele¬ 
scope/Very Large Telescope (NTT/VLT). The simulated 
orbits of S2 star are obtained by numerical integration 
of equations of motion where the hybrid gravitational 
potential is adopted, i.e. 

f = V, - - Vd) (r) , (27) 


where fi is the reduced mass in the two-body prob¬ 
lem. In that way we obtained the simulated orbit of 
S2 star around Galactic Centre in the weak field ap¬ 
proximation of hybrid gravity where eqs. (l24l) and (l25T l 
stand. Taking into account that y = y(0o, m^), the con¬ 
sidered weak field solution depends on the following 
three parameters: M, 0o, and m^. Mass M of the cen¬ 
tral object can be obtained independently using differ¬ 
ent observational techniques, such as e.g. virial analysis 
of the ionized gas in the central parsec (yielding 
M = 3 X lO^Msun), M - cr (mass - bulge velocity dis¬ 
persion) relationship for the Milky Way 1671] (yielding 
M - 9.4 X lO^Msun) or from Keplerian orbits of S-stars 
I 55n (yielding M = 4.3 x 10®Mj„„). Since our goal was 
not to make a new estimate of mass M using hybrid 
gravity, but instead to study the possible deviations from 
Keplerian orbit of S2 star which could indicate signa¬ 
tures for hybrid gravity on these scales, we adopted the 
last of three previously mentioned estimates for mass of 
the central object (M = 4.3 x lO^M^n), as well as the 
distance to the S2 star given by if^ (t/* = 8.3 kpc), and 
constrained only the remaining two free parameters ( 0 o, 
m^). Parameter 0o is dimensionless, while is given 
in AU“' (AU being astronomical unit), so that rep¬ 
resents a scaling parameter for gravity interaction. Non¬ 
zero values of these two parameters, if obtained, would 
indicate a potential deviation from GR. 
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In order to obtain the constraints on 0o and m^, these 
two parameters were varied. For each their combination 
the simulated coordinates x and y and velocity compo¬ 
nents Vx and Vy of S2 star were calculated. Calculations 
were performed for each observational epoch and then 
compared with its corresponding observed positions and 
velocities, between the observed and calculated co¬ 
ordinates of S2 star is minimized using LMDIFl routine 
from MINPACK-1 Fortran 77 library which solves the 
nonlinear least squares problems by a modihcation of 
Marquardt-Levenberg al gori thm 16411 (for more details 
on fitting procedures see 13611 1. 

4. Results: simulations vs observations 

Let us now discuss the numerical simulations that we 
want to compare with observations in order to select the 
range of the potential parameters (3.4). As we will see, 
analysis by hybrid gravity fixes better the observational 
data than the standard Keplerian analysis. 

4.1. Numerical calculation of S2 star orbit and orbital 
precession 

The simulated orbits of S2 star around the central ob¬ 
ject in hybrid gravity (blue solid line) and in Newtonian 
gravity (red dashed line) for 0o = -0.00033 and - 
-0.0028 (left panel), as well as for (pQ - -0.000033 and 
- -0.00028 (right panel) during 5 orbital periods, 
are presented in Fig. [1] As it can be seen from this fig¬ 
ure, hybrid gravity causes the orbital precession in the 
same direction as GR, but precession angle is much big¬ 
ger. When both (pQ and are decreased for an order of 
magnitude, the precession is much smaller (see the right 
panel of Fig. 1). This analysis also shows that Keplerian 
orbit is recovered when (po and tend to 0. 

We calculate orbital precession in hybrid modified 
gravity potential and results are reported in Fig. |2] as 
a function of 0o and m^. Assuming that the hybrid po¬ 
tential does not differ significantly from Newtonian po¬ 
tential, we derive the perturbed potential as 

F(r) = <l)(r)-0^(r) ; (I)^(r) = -— . (28) 

r 

The obtained perturbing potential is of the form: 

V(r) = [1 + (1/3) McPo/r. (29) 

1 + pQ 

and it can be used for calculating the precession angle 
according to Eq. (30) in Ref. IbSL 


-2L r z- dz dV (z) 

GMe^ J yi _ ff dz 


( 30 ) 


where r is related to z via: r - -. By differenti- 

l + ez 

ating the perturbing potential V(z) and substituting its 
derivative and (L = a(l - e^)) in the above Eq. (l30l l. 
and taking the same values for orbital elements of S2 
star like in Ref. lE^ we obtain numerically, for po - 
-0.00033 and - -0.0028, that the precession per or¬ 
bital period is 3°.26. 

Graphical representation of precession per orbital pe¬ 
riod for po in the range [-0.0009, -0.0002] and 
in [-0.0034, -0.0025] is given in the left panel of 
Eig. |2] As one can see, the pericenter advance (like 
in GR) is obtained. The precession per orbital pe¬ 
riod for pQ in the range [-0.0004, -0.0002] and in 
[-0.0029, -0.0027] is given in the right panel of Eig. |2] 


4.2. Comparison of theoretical results and observa¬ 
tions 

Let us give some constraints on parameters pQ and 
of hybrid gravity potentials according to current avail¬ 
able observations of S2 star orbit. However we should 
note that the present astrometric limit is still not suffi¬ 
cient to definitely conhrm that the S2 orbit deviates from 
the Keplerian one, but there is great probability that it is 
the case because the astrometric accuracy is constantly 
improving from around 10 mas during the first part of 
the observational period, currently reaching less than 
0.3 mas (see reference 1601]). There are also some re¬ 
cent studies that provide more and more evidence that 
the orbit of S2 star is not closing (see e.g. Eig. 2 in paper 
if^ L In this paper we fitted the NTT/VLT astrometric 
observations of S2 star, which contain a possible indica¬ 
tion for orbital precession around the massive compact 
object at Galactic Centre, in order to constrain the pa¬ 
rameters of hybrid gravity potential, since this kind of 
potential has not been tested at these scales yet. We 
have to stress that in the reference ifssll on page 1092, 
Eig. 13, authors presented the Keplerian orbit, but in or¬ 
der to obtain it they had to move the position of central 
point mass. In that way they implicitly assumed orbital 
precession. In our orbit calculation we do not need to 
move central point mass in order to get a satisfactory fit. 

In fact, our comparison with astronomical observa¬ 
tions represents upper bounds for precession angle on 
deviation from GR. The most probably results for pre¬ 
cession are in between these upper bounds and GR re¬ 
sults. In future, using more precise astronomical obser¬ 
vations, we could obtain more accurate results. 
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Figure 1: (Color online) Compaiisons between the orbit of S2 star in Newtonian gravity (red dashed line) and hybrid gravity during 5 orbital 
periods (blue solid line) for (left panel) = -0.00033 and = -0.0028, and for (right panel) = -0.000033 and = -0.00028. 
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Figure 2: (Color online) The precession per orbital period for 0o in the range [-0.0009, -0.0002] and in [-0.0034, -0.0025] (left panel), and 
00 in the range [-0.0004, -0.0002] and in [-0.0029, —0.0027] (right panel) in the case of hybrid modified gravity potential. With a decreasing 
value of angle of precession colors are darker. 
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Figure 3: (Color online) The maps of the reduced over the parameter space for all simulated orbits of S2 star which give at least the 

same or better fits than the Keplerian orbits. With a decreasing value of;^^^ (better fit) colors ingrey scale are darker. A few contours are presented 
for specific values of reduced given in the figure’s legend. 



Figure 4: (Color online) The same as in Fig. [3] but for the zoomed range of parameters. 
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Figs. |3]and|4]present the maps of the reducedover 
the 00 - parameter space for all simulated orbits of 
S2 star which give at least the same or better fits than the 
Keplerian orbits. These maps are obtained by the same 
fitting procedure as before. As it can be seen from Figs. 
|3] and |4] the most probable value for the parameter 0o 
in the case of NTT/VLT observations of S2 star is be¬ 
tween -0.0009 and -0.0002 and for the parameter is 
between -0.0034 and -0.0025 (see the darkest regions in 
Figs. IHl. In other words, we obtain reliable constraints 
on the parameters 0o and of hybrid modified gravity. 
The absolute minimum of the reduced - 1.503) 

is obtained for 0o = -0.00033 and - -0.0028, respec¬ 
tively. 

We simulated orbits of S2 star around the central ob¬ 
ject considering both the hybrid gravitational potential 
and the Newtonian potential. Our analysis shows that 
the hybrid modified gravity potential induces the pre¬ 
cession of S2 star orbit in the same direction of GR. We 
used these simulated orbits to fit the observed orbits of 
S2 star. The best ht (according to NTT/VLT data) is 
obtained for the 00 from between -0.0009 and -0.0002, 
and for the between -0.0034 and -0.0025. This range 
corresponds to scale parameter from (1/0.0034) AU 
to (1/0.0025) AU (=s 300 - 400 AU, i.e 1.4 - 1.9 mpc) 
which is comparable to the size of S2 star orbit. 

We believe that comparison with astronomical ob¬ 
servation is important, and data we used are the best 
currently published and available. GR predicts that the 
pericenter of S2 star should advance by 0°. 18 per orbital 
revolution 11551] . Using our fitting procedure, we get a 
much bigger precession 3°.26. Figure 2 in this paper 
gives theoretically calculated precession per orbital pe¬ 
riod for hybrid gravity of (po-m^ parameter space. In the 
future, with much more precise data maybe observation 
will hnd smaller value of precession, and using Figure 
2 (i.e. the same procedure), we will be able to get again 
hybrid gravitational parameters and 0o, hoping that 
observations will give smaller values. We calculated the 
map of parameters theoretically for broad range of pre¬ 
cession angles. More precise observations probably will 
change best fit parameters, but procedure for theoretical 
calculation will be the same. 


5. Conclusions 

In this paper, the orbit of S2 star around the galactic 
Centre has been investigated in the framework of the hy¬ 
brid modified gravity. Using the observed positions of 
S2 star, we constrained the parameters of hybrid modi¬ 
fied gravity. Our simulation results are; 


1 . the range of values for 0o parameter, coming from 
S2 star, is between -0.0009 and -0.0002; 

2. the range of is between -0.0034 and -0.0025; 

3. precession of S2 star orbit, in the hybrid modihed 
gravity potential, has the same direction as in GR, 
but the upper limit in magnitude is much bigger 
than GR. 


The above results allow to compare the orbital mo¬ 
tion of S2 star in the framework of hybrid gravity with 
analogous results in other theories. In particular, hy¬ 
brid gravity can be compared with metric f{R) mod¬ 
els, discussed in il29[ 13611 and with f{R, 0), discussed 
in 0. Also in these papers, the motion of S2 star has 
been studied according to the effective gravitational po¬ 
tentials achieved in the weak held limit. As discussed 
above, the main reason to introduce hybrid gravity lies 
on the fact that models like f{R) gravity (both in metric 
and Palatini formalism) and f{R, 0) gravity suffer prob- 
lems in passing the standard Solar S yste m tests ll4d.l4in . 
On the other hand, as reported in 15111 . hybrid gravity 
allows to bypass shortcomings deriving from local tests 
and connect models to galactic dynamics and late time 
cosmic acceleration. Using S2 star orbits, it is possi¬ 
ble to achieve additional constraints at sub-parsec scales 
and promote this model with respect to other extended 
gravity approaches. 

In particular, 0o and are the specihc parameters of 
hybrid gravity and differ from f{R) gravity models both 
in metric and Palatini formalism. In the case of f{R, 0) 
gravity, it is possible to achieve a Sanders like poten¬ 
tial; the parameter is also present and could have the 
same value, but the parameter 0o of hybrid gravity and 
a of f{R, 0) differ 1371]. The two effective gravitational 
potentials, in the week held limit, have similar, but not 
the same forms at sub-parsec scales. 

In conclusion, the comparison of the observed orbits 
of S2 star and theoretical calculations performed by the 
hybrid modihed gravity model can provide a power¬ 
ful method for the observational test of the theory, and 
for observationally discriminating among the different 
modihed gravity models. It seems that hybrid gravity 
potential is sufficient in addressing the problem of dark 
matter at galactic scales im, and it gives indications 
that alternative theories of gravity could be viable in de¬ 
scribing galactic dynamics. 

Furthermore, orbital solutions derived from such a 
potential are in good agreement with the reduced 
deduced for Keplerian orbits. This fact allows to hx 
the range of variation of 0o and m^. The precession 
of S2 star orbit, obtained for the best ht parameter val¬ 
ues (00 from -0.0009 to -0.0002 and from -0.0034 
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to -0.0025), has the positive direction, as in GR, but for 
these values of parameters, we obtain much larger or¬ 
bital precession of S2 star in hybrid gravity compared 
to prediction of GR. 


We can conclude that hybrid gravity effective po¬ 
tential is probably the best candidate among the other 
considered gravity models such as e.g. R" ||^, 381. 
Yukawa-like and Sanders-like to explain grav¬ 
itational phenomena at different astronomical scales. 


It is important to stress that our comparison with as¬ 
tronomical observations represents only upper bounds 
for precession angle on the deviation from GR. Al¬ 
though observational data seem to indicate that the S2 
star orbit is not Keplerian, the nowadays astrometric 
limits are not sufficient to unambiguously confirm such 
a claim. We hope that forthcoming observational data 
will allow more accurate measurements of stellar posi¬ 
tions. 


A final remark is due now. From an astrophysical 
point of view, the main motivation to introduce hybrid 
gravity is to address the problems of dark matter and 
dark energy EH. First of all, we have to say, accord¬ 
ing to the observations, that dark matter has very neg¬ 
ligible effects around the Galactic Centre . Despite 
of this fact, here we adopted hybrid gravity dynamics 
only to fit the orbit of S2 star around the Galactic Cen¬ 
tre. The interest of the reported results, if confirmed, 
lies on the fact that hybrid dynamics is independent of 
the dark issues but can be connected to a fine analy¬ 
sis of geodesic structure. In other words, the further 
gravitational degrees of freedom, coming from hybrid 
gravity, contribute to dynamics as soon as orbital anal¬ 
ysis related to GR is not sufficient to describe in detail 
peculiar situations as those around the Galactic Centre. 
However, in order to better confirm this statement, one 
needs more precise astronomical data describing stellar 
dynamics around Galactic Center. 
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